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ABSTRACT 


vA  model  for  the  prediction  of  the  elastic  response  of  composite  materials 
over  wide  ranges  of  concentration  Is  presented.  The  method  is  based  on  the 
mathematical  analogy  between  the  motion  of  particles  suspended  In  viscous  media 
and  the  elastic  deformation  of  solids.  The  system  Is  locally  described  by 
linearized  forms  of  the  elastic  moduli.  The  validity  of  these  relations  Is  then 
extended  to  all  concentrations  under  the  assumption  that  any  new  portion  of 
filler  "sees"  the  existing  composite  as  a  non-interacting  homogeneous  matrix. 

The  method  predicts  the  behavior  of  composites  with  solid  spherical  Inclusions 
and  foams  over  wide  ranges  of  concentration.  The  model,  free  from  adjustable 
parameters,  shows  excellent  agreement  with  existing  experimental  data.  The 
extension  of  the  method  to  other  Inclusion  geometries  Is  straightforward. 


By- - 7" 

Dlatrlt-utlci-'/ — ----- 

~AvyllaHltty  Codes^ 

-  ;Av^l  sr.d/or 


Special- 


INTRODUCTION 


In  the  search  of  a  mathematical  model  for  the  mechanical  behavior  of  hetero- 
geneous  systems,  a  complete  analogy  is  found  between  the  basic  equations  for 
elasticity  theory  and  those  for  the  motion  of  suspended  particles  In  viscous 
media.  The  analogy  can  be  easily  proved  by  writing  the  equations  of  motion  for 
an  incompressible  material  in  terms  of  displacements: 


(1) 


where  G  is  the  shear  modulus,  u^  is  the  displacement,  F^  Is  the  body  force,  p 
is  the  reactive  pressure,  and  p  Is  the  density. 

For  comparison,  the  Navier-Stokes  equations  for  an  Incompressible  Newtonian 
fluid  are: 


■p.t  *  '"'i *»Fi  •»  l  fj1-*  Vi,j)  <2> 

where  n  is  the  viscosity,  v  the  velocity,  and  F  the  body  force  per  unit  mass. 
Equations  1  and  2  have  a  term  by  term  equivalence,  with  the  exception  of  the 
nonlinear  component  v.v.  representing  a  convective  acceleration.  Under 
creeping  flow  conditions,  acceleration  terms  become  irrelevant  and  a  complete 
analogy  can  then  be  established.  Given  a  certain  geometty,  the  theoretical  equations 
for  a  viscous  flow 'problem  and  those  for  the  elastic  deformation  of  a  solid  are 
of  the  same  form: 

shear  rate  :  viscosity 
shear  strain  :  shear  modulus 


For  a  suspension  of  rigid  particles  In  an  Incompressible  matrix,  there  is  a 
simple  relationship  between  relative  viscosities  and  relative  shear  moduli 


(relative  to  the  unfilled  matrix): 


Thus,  If  a  theory  for  the  viscosity  of  a  filled  system  or  suspension  Is 
available.  It  can  be  used  to  estimate  the  shear  modulus  of  an  elastic  system. 

A  second  analogy  can  be  also  found  In  relation  to  the  methods  used  In  the 
solution  of  both  problems,  as  Indicated  In  Figures  1  and  2.  In  very  general 
terms.  It  Is  found  that  there  Is  an  Important  group  of  contributions  based  on 
the  establishment  of  an  Idealized  geometry  and  packing  arrangement  (top  of 
Figures  1  and  2).  A  second  group  Is  based  on  the  determination  of  bounds,  which 
avoids  (or  reduces)  the  Idealization  problem  by  use  of  energy  and  variational 
techniques  (bottom  of  Figures  1  and  2).  Typical  limitations  of  geometric  models 
and  bound  methods  gave  origin  to  a  third  group  of  contributions  containing  ele¬ 
ments  of  both  techniques  (middle  of  Figures  1  and  2). 

With  regard  to  the  problem  of  predicting  the  elastic  performance  of  com¬ 
posite  materials,  (Figure  2),  most  of  the  models  deal  with  specific  idealized 
geometries  and  packing  arrangements.  These  theories  provide  relations  between 
any  two  Independent  elastic  constants  among  a  shear  modulus  (G),  a  bulk  modulus 
(K)  and  a  Young's  modulus  (E): 

*  ■  W."  *«•  V  c> 

s  -G(GV  K,.  G,,  K,.  c) 

where  subindices  m  and  i  refer  to  the  matrix  and  the  Inclusion,  and  c  Is  the 
volume  fraction  of  the  Inclusion.  Smith  [11]  noted  that  through  algebraic 


manipulation,  most  of  these  theories  can  be  displayed  In  the  form  of  a  genera* 
11  zed  "rule  of  mixtures": 


G  ’  Gm  +  (G1  “  V*  "  (l  "  ♦)<*m  +  *Gi 

where  #  takes  different  functional  forms.  According  to  the  most  familiar 
contributions, 

t  _  (8  -  10,J8.  -  S.J8.  _ 

(8  -  lOvJG,  *  (7  -  5,B)Gm 

(approximate  Van  der  Poel)  [11,16] 

(8  -  IOvJG  +  (7  -  5v_)G 
$  a  _  ffl  ffl  in  £ 

(8  -  10vm)G,  *  (7  -  5vb)Gb 

(Kerner)[17] 

(8  -  10v)G  +  (7  -  5v )G  c 
*  "IT-  10v)S1  ♦  (7  -  5v)G 

(Bud1ansky)[18] 

where  G  and  v  are  the  shear  modulus  and  Poisson's  ratio  of  the  composite. 

All  of  the  previous  models  developed  by  elastic  analysis,  l.e.,  deformation 
of  the  matrix  and  inclusion,  have  a  common  feature:  a  specific  Idealized 
geometry  and  packing  arrangement  has  to  be  assumed.  In  many  composite 
materials,  this  regularity  In  geometry  does  not  exist,  and  the  material  cannot 
be  Identified  with  any  particular  arrangement  or  ensemble  valid  throughout  the 
whole  structure.  Thus,  the  previous  treatments.  In  the  strictest  sense,  do  not 
truly  represent  a  real  composite.  One  way  to  avoid  this  Idealization  problem  is 


through  the  use  of  variational  techniques,  which  do  not  need  to  use  a  precise 
geometric  model. 

The  most  elementary  type  of  bounds  can  be  calculated  through  the  formalism 
of  the  complementary  energy  and  strain  energy  theorems  (Voigt  [12].  Paul  [13]). 
However,  the  resulting  bounds  are  often  too  far  apart  to  be  useful. 

In  order  to  improve  these  bounds,  Hashin  and  Shtrikman  [14]  formulated  new 
variational  theorems  on  the  basis  of  polarization  tensors.  However,  limita¬ 
tions  In  the  applicability  of  bounds  still  persist.  In  general.  It  Is  not  known 
if  the  calculated  bounds  are  Indeed  the  most  restrictive,  and  the  problem 
becomes  a  trade-off  between  the  amount  of  Information  Introduced  In  the 
formulation,  and  the  resulting  contiguity  of  the  bounds. 

As  a  rule.  It  may  be  said  that  the  usefulness  of  bounds  Is  restricted  to 
systems  where  the  properties  of  the  phases  are  similar,  i.e.,  metal  alloys.  In 
polymer  reinforcement,  however,  the  differences  between  matrix  and  Inclusion  may 
be  considerable.  As  a  consequence,  the  upper  and  lower  bounds  tend  to  separate, 
making  this  technique  Inadequate  for  any  predictive  purpose. 

PROCEDURE 

As  Indicated  earlier,  a  basic  problem  In  the  theory  of  viscosity  of  suspen¬ 
sions  Is  the  determination  of  correlations  useful  at  high  concentrations.  Most 
of  the  contributions  are  based  oh  extensions  of  Einstein's  law  of  viscosity  at 
Infinite  dilution. 

The  analogous  problem  of  correlating  the  elastic  properties  of  heterogeneous 
systems  over  wide  ranges  of  concentration  has  not  been  successfully  solved.  As 


Mill  be  discussed  later,  predictions  according  to  different  theories  shoM  a  con¬ 
siderable  lack  of  agreement. 

The  proposed  model  Is  based  on  a  differential  scheme  similar  to  that  used 
by  Brinkman  [3]  and  Roscoe  [2]  for  the  prediction  of  the  viscosity  of  suspen¬ 
sions  at  high  concentrations.  In  a  later  publication,  Roscoe  [9]  tried  to 
extend  the  method  to  the  prediction  of  elastic  properties  of  composites. 
McLaughlin  [30]  presented  a  similar  differential  scheme  and  evaluated  Its  esti¬ 
mates  In  relation  to  the  Hashln-Shtrlkman  bounds. 

The  present  model  is  based  on  the  calculation  by  elastic  analysis  of  any 
two  independent  parameters  for  the  composite,  namely  the  shear  modulus,  the  bulk 
modulus,  or  the  Young's  modulus.  Calculations  may  proceed  as  follows  [29]: 

1)  describe  the  deformation  state  (simple  shear,  etc.) 

2)  choose  adequate  coordinates  for  a  given  geometry, 

3)  state  the  deformation  equations  for  matrix  and  Inclusion, 

4)  state  continuity  of  stress  at  the  Interface, 

5)  reduce  volume  Integrals  to  surface  Integrals  by  Eshelby's 
formula  [5], 

6)  Integrate. 

For  spherical  Inclusions  and  dilute  conditions,  the  shear  modulus  results  as: 


and  the  bulk  modulus  Is: 


(1/Km)  (*1 -Km) 


1  +  gc 


(4) 


It.  1  ♦  cot,  -  «,)/<«.  ♦  (4/3JG.J 
where  c  notes  the  volume  fraction,  and  the  subindices  m  and  1  refer  to  the 
matrix  and  the  Inclusion. 

For  an  Infinitely  rigid  Inclusion  (G^  •)  in  an  Incompressible  matrix 

(v^  ■  0.5),  equation  (3)  results  in  a  relation  similar  to  Einstein's  law  of 
viscosity: 


jL-  l  +  2.5  c  (5) 

nm  Gm 

Consider  a  suspension  containing  a  small  concentration  (volume)  Cj  of  spheres 
into  which  Is  placed  a  small  concentration  c2  of  spheres,  and  Into  this  again 
a  small  concentration  c-j  of  spheres,  and  so  on,  up  to  n  sets  of  spheres,  and  let 
the  total  concentration  be  c.  Then  on  adding  another  set  of  spheres  of 
concentration  (volume  fraction)  cn+^,  the  new  total  concentration  (volume 
fraction)  Is: 

C(1  "  cn+l}  +  cn+l  (fi) 

ft 

So  the  Increase  In  concentration  Is: 


“  ■  'final  '  'Initial  *  W1  *  c> 


(7) 


If  the  original  suspension  Is  considered  a  homogeneous  medium  (matrix)  around 
the  new  spheres  which  are  In  small  concentrations.  It  Is  possible  to  express 


the  Increment  In  shear  modulus  according  to  equation  (3)  as: 


Defining 


^  ■ 1  *  fVi 

n 


AG  *  Gn+1  “  Gn 


and  replacing  cn+1  according  to  equation  (7), 


M  »  f  c  ,  - 
G  «+l  1-c 


For  Infinitely  small  Increments  In  concentration. 


ac  ♦  dc 
aG  ♦  dG 


«  -Gf  (9) 

dc  (1-c)  ’ 

This  differential  expression  for  the  Increment  In  concentration  can  now  be 
used  in  equation  (3).  After  following  a  similar  reasoning  for  the  bulk  modulus 
K,  equations  (3)  and  (4)  become: 


dG  _  *  G  lSU-vKl-G^/G] 
dc  *  [7-5v  +  2(4-5v)G1/G](l-c) 


(10) 


a _ (Ks  -  K) _ 

dc  (1  ♦  C(Ki_K)/(K+(4/3)G) (1-c) 


where  v  »  iiJLl_2£L.  i s  the  Poisson's  ratio. 

2(3K  +  G) 

Equations  (10)  and  (11)  constitute  a  coupled  system  to  be  numerically 
Integrated  with  the  following  boundary  conditions: 


c  -  0 

G 

“Gm 

(matrix  properties) 

K 

aKm 

c  -  1 

G 

*G1 

(inclusion  properties) 

K 

*K1 

* 

RESULTS 

1)  Solid  Inclusions 

Smith  [11]  compared  the  behavior  of  well  known  models  for  the  prediction  of 
the  mechanical  response  of  composites,  i.e..  Van  der  Poel  [16],  Kerner  [17], 
Budlansky  [18],  etc.  In  order  to  test  the  performance  of  the  present  model,  the 
same  parameter  values  were  used  In  the  simulation.  Data  In  Figure  3a  have  been 
calculated  with  a  shear  modulus  ratio  G,/G_  ■  30  and  Poisson's  ratios  v*  *  0.25 

i  m  1 

and  vm  *  °*40«  representing  the  case  of  glass  spheres  embedded  in  an  epo*y 
matrix. 

Figure  3b  has  been  calculated  with  G./Gm  ■  70,000,  v.  *  0.25,  v  -  0.5, 

•  m  1  m 

representing  the  expected  shear  properties  of  glass  spheres  Imbedded  in  a 
lightly  vulcanized  matrix  of  natural  rubber. 


In  both  cases,  a  considerable  discrepancy  among  the  different  models  is 
observed.  The  inconveniency  of  bounds  calculated  according  to  Hashin  and 
Shtrikman  [14]  is  also  stressed.  The  predictions  of  the  present  simulation  are  close 
to  those  of  Van  der  Poel  [16]. 

Richard  (7)  presented  experimental  data  for  elastic  modulus  and  Poisson's 
ratio  of  a  polyester  matrix  reinforced  with  glass  spheres  (G^  »  5.94  108  Pa, 
vffl  *  0.45,  G^  *  2.90  1010  Pa,  *  0.21).  This  author  compared  his  experiments 
with  the  predictions  of  several  theories.  Figure  4  Indicates  the  excellent  pre¬ 
diction  of  the  present  model,  while  theories  by  Kerner  [17]  and  Hill  [20,21] 
only  perform  well  at  filler  concentrations  below  20  percent  volume. 

Richard's  data  [7]  were  also  used  by  Smith  [11]  for  comparison  with  other 

theoretical  contributions  as  indicated  in  Figure  5a  for  relative  moduli,  and 

Figure  5b  for  Poisson's  ratio.  Calculations  were  performed  with  Gm  *  5.84 

108  Pa,  »  0.44,  G.  »  2.90  101°  and  v  »  0.21. 
m  i  1 

Smith  [23]  compared  his  own  experimental  findings  on  epoxy  resin-glass 
spheres  systems  with  several  theoretical  predictions.  The  properties  of  matrix 
and  inclusion  are  E^  *  2.68  109  Pa,  vm  *  0.394,  Ei  s  7.6  1010  Pa,  vi  *  0.23. 

Results  from  those  predictions  are  plotted  and  compared  with  the  present  model 
in  Figure  6. 

2)  Foams 

Equations  (10)  and  (11)  have  been  Integrated  for  the  case  of  a  foam  by 
assuming  negligible  moduli  for  the  Inclusion, 


The  simulation  was  compared  with  the  experimental  results  of  Gent  and  Thomas 
[24]  for  a  vulcanized  rubber  foam  (Em  ■  2.59  10®  Pa,  \>m  *  0.49).  Figure  7  shows 
excellent  agreement  in  the  prediction  of  relative  modulus  versus  relative  den¬ 
sity  6, 

4  =  ~6..foam  ,  l  -  c  (13) 

6  rubber 


where  c  is,  as  before,  the  volume  fraction  of  the  inclusion. 

The  performance  of  the  model  in  the  prediction  of  foam  properties  was  also 
tested  with  the  experimental  work  of  Moore  et  al.  [25]  for  a  variety  of  foamed 
thermoplastics  (PVC,  polypropylene  copolymer,  styrene-acrylonitrile  copolymer, 
etc.).  These  authors  observed  that  plots  of  relative  modulus  (relative  to  the 
solid  matrix)  versus  relative  densities  (6)  could  be  correlated  with  a  single 
curve  in  shear,  tension  and  compression  according  to: 


'foam 

:solid 


[: 


density  of  foam 
density  of  solid 


l 


(«)2  »  (1  -  C)2 


(14) 


Figure  8  compares  the  predictions  of  several  theories  with  the  empirical 
findings  of  Moore  et  al.  [25]  and  the  results  of  the  present  model.  Simulations 
were  performed  for  a  polypropylene  copolymer  characterized  by  Em  *  1.13  GPa, 
vm  *  0.41.  Again,  the  proposed  method  shows  total  agreement  with  the  experimen¬ 
tal  values. 

3)  Low  modulus  Incompressible  inclusions  (rubbers) 

For  the  limiting  case  of  an  Incompressible  matrix  (v^  ■  0.5)  with  G^  *  0, 


equation  (10)  admits  an  analytical  solution  as 

p-  *  (1  -  c)5/3  (15) 

Gm 

which  represents  the  behavior  of  rubber  modified  polymers. 

Holliday  et  al.  [31]  compared  several  relevant  theories  for  a  case  repre¬ 
senting  soft  spheres  In  a  rigid  matrix  (Figure  9).  The  experimental  data 
correspond  to  acrylonltrl le-butadl ene-sty rene  (ABS)  polymers  which  consist  of  a 
dispersion  of  soft  polybutadiene  rubber  spheres  in  a  styrene-acrylonitrile  (SAN) 
copolymer  matrix  which  has  a  relatively  high  modulus. 

The  different  theories  examined  are  not  In  good  agreement.  The  proposed 
model  was  tested  with  typical  values  for  styrene-acrylonitrile  copolymer  (SAN) 
(Em  *  3.1  10" ^  MPa,  vm  ■  0.38)  [32],  and  /Em  *  0.  Results  presented  in 
Figure  9  Indicate  good  prediction  of  the  experiments. 

CONCLUSIONS 

Despite  many  contributions,  the  problem  of  predicting  the  elastic  behavior 
of  composite  materials  over  wide  ranges  of  concentration  has  not  been  satlsfac- 
torlly  solved.  .  . 

Different  theories  based  on  well  defined  geometric  models  present  a  con¬ 
siderable  divergence  of  results.  The  use  of  bound  methods  has  been  proven  Ina¬ 
dequate  for  problems  In  polymer  reinforcement,  where  solutions  become  a  trade¬ 
off  between  the  contiguity  of  the  bounds  and  the  amount  of  Information 
(complexity)  Introduced  In  the  model. 


In  the  present  work,  a  correct  linearized  relation  for  the  elastic  defor¬ 
mation  of  matrix  and  Inclusion  Is  assumed  at  each  concentration.  This  relation, 
which  Is  valid  locally,  can  then  be  extended  to  all  concentrations  under  the 
assumption  that  any  new  portion  of  filler  “sees"  the  existing  composite  as  a 
non-interacting  homogeneous  matrix. 

The  model,  presenting  no  adjustable  parameters,  accurately  represents 
available  experimental  data  on  solid  spherical  Inclusions  and  foams.  The  method 
can  be  readily  extended  to  other  geometries  for  the  Inclusion. 
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FIGURE  CAPTIONS 


Figure  1  Models  for  viscosity  of  suspensions. 

Figure  2  Models  for  prediction  of  elastic  response  In  composite  materials. 

Figure  3a  Relative  shear  modulus  versus  volume  fraction  of  filler  for  G^/G^  * 

30,  Vj  ■  0.25,  vm  ■  0.4  (from  Smith  [11]).  Solid  curves,  model 
performance;  dashed  curves,  predictions  with  different  models: 

(1),  approximate  Van  der  Poel  [16];  (2),  Kerner  or  Hashln  and 

Shtrlkman  highest  lower  bound  [14];  (3),  corrected  Van  der  Poel 

[16];  (4),  Budlansky  [18]. 

Figure  3b  Relative  shear  modulus  versus  volume  fraction  of  filler  for  G./G  » 

1  m 

70,000,  v,j  *  0.25,  m  *  0.5  from  Smith  [11],  Solid  curve,  model 
performance;  dashed  curves,  predictions  with  different  models. 

(1),  approximate  Van  der  Poel  [16];  (2),  Hashln  and  Shtrlkman 

highest  lower  bound  [14];  (3),  corrected  Van  der  Poel  [16]; 

(4),  Budlansky  [18], 

Figure  4  Elastic  modulus  versus  volume  fraction  of  filler  for  the  case  of 

glass  beads  In  a  polyester  matrix.  Experimental  data  from  Richard  [19]. 
Solid  curve,  model  performance;  dashed  curves,  predictions  according 
to  1),  Isostrain;  2),  Isostress;  3),  Kerner  [17];  4),  Hill  [20,21]. 


Figure  5  Relative  Young's  modulus  E/Em  (5a)  and  Poisson's  ratio  (5b)  versus 
volume  fraction  of  filler  for  the  case  of  glass  beads  In  a  polyester 
matrix.  Experimental  data  from  Richard  [19].  Solid  lines,  model 
performance;  dashed  lines,  predictions  with  different  models. 

1),  Hashln  and  Shtrlkman  highest  lower  bound  [14];  2),  Hashln  and 
Shtrlkman  least  upper  bound  [14];  3),  Van  der  Poel  [16];  4), 

Budlansky  [18]. 


Figure  6  Relative  Young's  modulus  E/Effl  versus  volume  fraction  of  filler  for 
the  case  of  glass  spheres  In  an  epoxl  matrix.  Experimental  data  by 
Smith  [23].  Solid  curve,  model  performance;  dashed  curves,  predic¬ 
tions  with  models  by  1),  Hashln  and  shtrlkman  highest  lower  bound  [14]; 
2).  Hashln  and  Shtrlkman  least  upper  bound  [14];  3),  Van  der  Poel 

[16];  4),  Budlansky  [18]. 

Figure  7  Young's  modulus  of  a  highly  vulcanized  rubber  foam  Ef  relative  to  the 
solid  rubber  Em  versus  volume  fraction  of  rubber  In  the  foam.  Solid 
line,  model  performance.  Experimental  points  were  determined  by  Gent 
et  al.  [24]. 

Figure  8  Relative  moduli  of  foams  versus  relative  densities.  Calculations 

here  have  been  performed  for  the  thermoplastics  described  by  Moore  et 
al.  [25].  Solid  curves,  model  performance  for  1),  Young's  modulus; 

2) ,  shear  modulus.  Dashed  curves,  predictions  according  to  different 

models.  3),  Kerner  [17];  4),  square-1  n-square  model  [27];  5),  cube- 

in-cube  model  [27];  6),  law  of  mixtures.  Moore  et  al.  [25]  correlated 

experimental  moduli  with  a  square  law  which  results  coincident  with 
curve  1  (model  prediction). 

Figure  9  Relative  modulus  of  acrylonltrlle-butadlene-styrene  (ABS)  copolymer 
versus  volume  fraction  of  polybutadiene.  Experimental  data  by 
Holliday  et  al.  [31].  Solid  curve,  model  performance.  Dashed  curves, 
predictions  with  different  models.  1),  Van  der  Poel  [16];  2)  Parallel; 

3) ,  Kerner  [17];  4),  Relner-Hashln  [33];  5),  Series. 
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